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Abstract - We investigate the topological zeta function for unimodal maps in 
general and dynamical zeta functions for the tent map in particular. For the generic 
situation, when the kneading sequence is aperiodic, it is shown that the zeta functions 
have a natural boundary along its radius of convergence, beyond which the function 
lacks analytic continuation. We make a detailed study of the function rinLo(l — 
z 2 ) associated with sequences of period doublings. It is demonstrated that this 
function has a dense set of poles and zeros on the unit circle, exhibiting a rich 
number theoretical structure. 



1 Introduction and preliminaries 

In chaotic dynamics one aims at a global description of the phase space flow rather than a 
description of individual trajectories. The appropriate tool for this analysis is evolution operators 
(Ruelle-Perron-Frobenius operators). It requires a rather extensive mathematical apparatus to 
understand the spectrum of this operator in detail. The problem is that one has to restrict 
the functions space, acted on by the operator, in order to capture only the physically relevant 
spectrum. However it seems as if, in many cases, this relevant spectrum can be identified with 
the spectrum of poles of zeta functions or (even better) zeros of Fredholm determinants. These 
functions require only knowledge of the periodic orbits and their invariants, such as stability and 
length and no reference to specific function classes. The link between the spectra of evolution 
operator and their zeta functions has been rigorously established only for restricted classes of 
systems, like e.g. expanding (or hyperbolic) systems with a symbolic dynamics of finite subshift 
type [EL ^|. For such systems the Fredholm determinants may even be entire functions. 

However, there are many ergodic systems of physical relevance that do not fall into this class 
of systems. We can then no longer expect singularities to be absent from the zeta function. 
There are strong indications that for intermittent systems the zeta functions exhibit branch point 
singularities & 0, |j| . The corresponding branch cut could possibly be associated with a continous 
component of the spectrum of the evolution operator. There is evidence that these singularities 
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do carry important information about the dynamics of the system, such as power law decay of 
correlations and phase transitions among e.g. the generalized Lyapunov exponents j|, ||, 0]. 

In this paper we will study the analytic structure of the zeta function for expanding systems 
with no symbolic dynamics of finite subshift type. 

We will study the simplest conceivable system, namely one-dimensional unimodal (one- 
humped) maps x i— ► fx (x) = A • g(x) with one external control parameter A ||, || [Tc[ |. 

The kneading sequence £ A is the orbit of the critical point, represented by a string of R and 
L's depending on the branch of the map visited. The ordering of kneading sequences with respect 
to the parameter A do not depend on the details of the unimodal g(x) if it has strictly negative 
Schwartzian derivative |10, S|. For the tent map this derivative is zero and the sequences of period 



doublings arc squeezed into a zero length intervals in parameter space. 
The kneading sequence can be of one of three types 

a) It maps to the critical point again, after n iterations. If so, we adopt the convention to 
terminate the kneading sequence with a C, and refer to the kneading sequence as finite. 

b) Preperiodic, i.e. it is infinite but with a periodic tail. 

c) Aperiodic. 

Much dynamic information is encoded in the following one-parameter family of zeta functions 
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The product in (|l|) runs over all primitive periodic orbits p, having period n p and stability 
A p = d ^ d J \ x =x p with x p being any point along p. We will only study a very special case of this 
function in detail, namely the topological zeta function 

l/C top (z) S l/C(/3-l,2) = n(l-^ P ) > (2) 

p 

and then turn to some speculations about f3 < 1 at the end. The leading zero zq of l/Ctop(z) 
(the one with smallest modulus) and the topological entropy h are related by h = — logzo. There 
is one unimodal map for which l/Q op (z) contains all information concerning l/£(/3, z) for all 0, 
namely the following piecewise linear map, called the tent map [[l2| [l3| [l4| 

fix) - I X ' X [0,1/2] 

x~m- \ X -(i-x) xe (1/2,1] ■ { } 

where the parameter A £ [0, 2] Since |A P | = A" p we have 

l/C(P,z) = l/Cto P (z/\ 1 - ) ■ (4) 

The topological entropy is h = — log A. 

The set of periodic points of the tent map is countable. A consequence of this fact is that 
the set of parameter values for which the kneading sequence is finite or eventually periodic are 
countable and thus of measure zero and consequently the kneading sequence is aperiodic for almost 
all A. For general unimodal maps the corresponding statement is that the kneading sequence is 
aperiodic for almost all topological entropies. 

It is, in fact, easy to write down the expanded topological zeta function corresponding to a 
given kneading sequence J A = PC, of length n, where P = S1S2 ■ ■ ■ s n _i is a string of symbols 
where each symbol = L or R. Now let a L = 1 if Sj = L and aj = — 1 if Sj = R. In appendix A 
we show the following 
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Theorem 1 The expanded zeta function is a polynomial of degree n 

71-1 

i/Cto P (z) = l[(i-z;) = (i-z)-J2 b ^ i . ( 5 ) 

p i=0 

w/iere 

6o = 1 /g\ 

bi = bi-ia,i 1 < i < n — 1 

J/ i/ie kneading sequence is infinite (i.e. type b) or c)) the finite sum is replaced by an infinite 
series. 

The function YllZo biz 1 is called the characteristic polynomial in the literature [t5). In 
ref. fl2|| it appeared as a characteristic function whose zeros corresponds to eigenvalues of the 
evolution operator acting on the space of piecewise constant functions. 

An important consequence of ([|) is that the sequence {b{\ has a periodic tail if and only f 
the the sequence {c^} has one. We will make use of this observation in section 3. 



2 A quick journey through parameter space 

The universal sequence of finite kneading sequences is a well studied subject but to make ourselves 
acquainted with (||) we will discuss some special cases . The letters a) to d) refer to table 1. 

a) J_ x = H°°(R), i.e. the harmonic extension of RC, see ||. This is the accumulation point 
of the first cascade of period doublings. The topological zeta functions is given by l/Q op (z) = 
(1 — z) • IIriLo(l — 2,2 )■ The kneading sequence is aperiodic. We will study this function in some 
detail in the next section. Similar cascades of period doublings accompanies all periodic orbits, 
or equivalcntly phrased, a sequence of harmonics follows every finite kneading sequence. 

b) / A = RLR°° . This is a so called band merging point. The inverse topological zeta function 
is l/(to P {z) = (1 — z)(l — 2z 2 )/(l + z). The kneading sequence is of type b). The preperiodic 
kneading sequence give rise to the denominator (1+z) above and thus a pole of l/Q op (z). The two 
leading zeros z = ±l/v2 are of the same size. For the tent map, a slight increase of the control 
parameter will yield a gap between the leading and next to leading zero and strong mixing. This 
point, and higher band merging points, has been extensively studied in the literature [ jl2[ p"3[ fblj 

c) T_ x = RLC. The last Sarkovski orbit enters the system, orbits of arbitrary period are 
allowed. The kneading sequence is RLC and via eqs (||) and (||) we get l/C,t op {z) = (1 — z) ■ (1 — 
z — z 2 ). This can also be realized as follows (cf. refs. jl9| for details). The kneading sequence 
is RLC, so the only forbidden subsequence is RLL. All allowed periodic orbits, except L, can 
can be built from a alphabet with letters RL and R, following ref. we write this alphabet as 
\ RL , R; L}, yielding the topological zeta function l/Q op (z) = (1 — z)(l — z — z 2 ). The leading 
zero is the inverse golden mean Zq = (y/E — l)/2. 

d) Xa = RL°° . The symbolic dynamics is an unrestricted binary one with alphabet {R,L}, 
the zeta function is simply l/Q op (z) = 1 — 2z. The kneading sequence is of type b) but the 
associated denominator (1 — z) gets canceled by the prefactor in (|J). 



3 Lacunary zeta functions 

Consider a function f{z) = ^CiZ 1 defined by a power series. A meromorphic functions has a 
well defined analytic continuation in the entire complex plane, except at the poles. Functions 
exhibiting branch point singularities can be continued except on certain imposed cuts. 
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RLRRRLRC 




a)H°°(R) 


EUot 1 - Z ) 


RLRRRC 




RLRRRRRC 




b) RLR°° 


(l-2z 2 )/(l + z) 


RLRRRRRRC 




RLRRRRC 




RLRRC 




it i \ /i 




c) RLC 


(l-z-z 2 ) 


RLLRLC 




RLLRLRC 




RLLRC 




RLLRRRC 




RLLRRC 




RLLRRLC 




RLLC 




RLLLRLC 




RLLLRC 




RLLLRRC 




RLLLC 




RLLLLRC 




RLLLLC 




RLLLLLC 




d) RL°° 


(l-2z)/(l-z) 



Table l: Ordered kneading sequences up to length seven and some longer 
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But there is a class of functions that cannot be analytically continued outside a region (typ- 
ically a disk), the boundary of this region is called a natural boundary. There are some special 
cases where it is easy to detect the presence of a such a boundary. The classical example is 
gap series or lacunary series, first studied by Fredholm, Weierstrass and Hadamard. The easiest 
example is perhaps provided by the series f{z) — ^z 2 . The increasing gaps in the sequence of 
coefficients make the function noncontinuablc beyond the circle of convergence \z\ = 1. 

Another example is when the coefficients can take on only a finite number of values and the 
sequence of coefficients does not ultimately become periodic, then the function again has the unit 
circle as a natural boundary Jig] . This result can be directly applied to the expansion of 1/Ctop, 
cf eq. (^), for an aperiodic kneading sequence. We can thus get the following theorem for free 

Theorem 2 The topological zeta function 1/Ctop for unimodal maps has the unit circle as a 
natural boundary for almost all topological entropies and for the tent map for almost all X. 

Below we are going to study a specific example of zeta functions exhibiting a natural boundary. 
The concept of harmonic extension, and the associated cascades of period doublings is a central 
concept for the description of unimodal maps. This motivates a close study of the function 

oo 

so^nt 1 -^") ■ ( 7 ) 

n=0 

The expansion of begins as S(z) = 1 — z— z 2 +z 3 — z 4 + z 5 . . .. The radius of convergence is 
obviously unity. The simple rule governing the expansion will effectively prohibit any periodicity 
among the coefficients making the unit circle a natural boundary. We will not follow this line of 
thought, we will instead make a detailed study of the fine structure of this boundary and, as a 
by product, conclude that the unit circle is opaque. 

First we rewrite S as 

oo 

S(z) = n(l- X ») , (8) 

71=0 

where the x„'s are defined iteratively: 



Xn = xl-1 

Xo = z 



(9) 



We write the Xn on polar form Xn = r n exp(2irix n ) where r n and x n are real numbers and x n is 
restricted to the unit interval. The iteration rule is 

2 2 n 

i n = 2x n -i mod 1 

Suppose now we choose rg close to, but inside the unit circle ro = 1 — e. Then r n will be close to 
unity for n < N ps log 2/ log e and then decrease to zero quickly. By choosing e we effectively select 
N(e) factors of the product representation of 2 and by decreasing e we can conclude whether 5(z) 
diverges to oo or or perhaps converge to a finite value as z approach exp(27ria;) from within; 
we thus refer to z = exp(27rix) as a pole, zero or a regular point respectively. This reasoning is 
made rigorous in appendix B. 

We will now focus our attention on the second recursion relation in ( |l7j| ) which happens to 
be the binary shift map, and we will restrict our attention to rational xq. Take xq — p/q to be 
proper reduced fractions, that is q £ N and 



p E fl q = {p; p e N, 1 < p < q, 0, q) = 1} 



(11) 
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oo 
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oo 
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oo 
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oo 
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2 ,l 2 ,oo 2 
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11 


20 
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0,oo 



Table 2: Some number theoretic data used to study the structure of the unit circle for 
the function Explanations given in section 3. 



where (p, q) denote the greatest common divisor of q and p. The number of elements in fl q is 
Euler's function <p(q). We are thus led to study the following map from fl q to itself 

p n = 2p n -i mod q . (12) 

We will separately consider the three cases, 
i/ q = 2 m 

ii/ q = 2"Y where (q J ,2) = 1 
iii/ (q, 2) = 1 

We begin with i/. If we set r = 1~ then only the first m factors of the product (Q) are nonzero 
and all the subsequent will be zero and the corresponding point on the unit circle z = exp(2irix) 
is a zero. These are in a sense trivial zeros and we will eventually discover that they are not the 
only ones. 

Case ii/ reduces to case i/ after m applications of the map. 

Now to the more interesting behaviour of the map ( |l2] ) for case iii/. Below we are going to 
apply some fundamental results from number theory, they can be found in any standard textbook. 
The text below is rather compact so in order to digest it could be advisable to verify table 2, step 
by step during the reading. 

Now 2 G Q q . It is well known that the set £l q is an abelian group with respect to multiplication 
modulo q. This has some important consequences. The map will move points in Q q along cycles, 
because the existence of transients would imply nonexistence of unique inverses. Secondly if the 
set fig is divided into several subcycles, all of them has the same length l q , and l q is a divisor of 
</>(<?). The number l q will be the smallest solution to 

2 l " = 1 mod q . (13) 

Due to Euler's theorem we see that l q = <f>{q) is always a solution but it is not necessarily the 
smallest. In number theory one says that 2 is a primitive congruence root modulo q if l q = 4>{q), 
or equivalently, that 2 is a generator of the residue class modulo q. 



Effect of pruning 



7 



First, if q — p is a prime then 2 may or may not be a generator. There is no other way 
than to check it by hand (or look it up in a table). For p < 100 we can use 2 as a generator for 
p = 3,5,11,13,19,29,37,53,59,61,67 and 83 but not for 7,17,23,31,41,43,47,71,73,79,89 and 97. If 2 
generates the residue class for q = p it also acts as generator for q = p a unless 2 P ~ 1 = 1 mod p 2 . 
This is an accident which is very rare if possible, it does not occur for any prime p < 1000. When 
building up table 2 we do not have to worry about this complication. 

Next suppose q is composite. Than, there does not exist any primitive congruence root and 
and l q < <f>(q). Indeed, if q = q\ ■ q<i and = 1 (and q%, q% ^ 1) it is fairly easy to show 

that l q < 4>{q)/ (4>(qi), ^(92)) < </ ) ('?)/2. We can be more explicit than so: decompose q into prime 
factors q = p^p^ 2 ■ ■ -Pr r where all on > 0. If all residue classes modulo p^ can be generated by 
2 then one can show that l q — {(^(p" 1 ), </>(.P2 2 ), ■ ■ ■ ,^(p? r )}, where {a, b} stands for least common 
multiple with obvious generalization for several arguments. 

The reader can now verify the third column of table 2. 

The analysis above can tell whether or not the set fi 9 splits up into subcycles under action 
of the map p > 2p mod q. We can also tell how many subcycles in total having period n. The 
reason is that the binary shift map admits a binary coding of its periodic orbits. There are 

M n = - V fi(n/d)2 d (14) 
n * — ' 

d\n 

primitive cycles of length n, /j is the Mobius function, see e.g. ref. jl9| for a derivation. For 
instance, Af 3 = 2 so there are two period three cycles, both corresponding to q = 7, cf. table 2. 
M4 = 3 so there are three cycles of period four, corresponding to q = 5, q = 9 and q = 15. 

Next we address the question if a given (sub)cycle c corresponds to a pole, a zero or a regular 
point of the function E(z). For cycle c we define 

^ c = Y[(l-e 2nip i^) . (15) 

From the previous reasoning it is clear that if \tp c \ < 1 then the points z — exp(27ri Pj/q c ) for all 
Pj G c are zeros of S(z). If ip > 1 they are poles and if ip = 1 they are regular points. It is not 
straight forward to express ip c (without numerical computation) for a given cycle but the product 
over all cycles with a given q c = q 

*,= n^= II (l-e 27 "^F g (l) (16) 

qc=q l<3<q,(j,q)=l 

is easily expressed. F q (x) are the so called cyclotomic polynomials frequently studied in number 
theory. We can again use well established results |H from number theory and state that 

If q is a prime power q = p a then ty q = p, and if q contains two or more distinct prime 
factors then ^ q — 1. 

The reader can now verify the fourth column of table 2. 

We were interested in the contribution ij) c from one cycle and if l q < 4>{q) we are only able 
to compute the product ^ q of contributions from all cycles corresponding to q. In some cases 
we find that two cycles c\ and C2 are related by complex conjugation (again a purely number 
theoretical property): = * C2 and thus |\P C1 | = l^cj. This is the case for e.g. q = 7. When 
there is no such symmetry there seems to preference to produce some c such that \^ c \ < 1 leading 
to zeros of S(z) that are surely non trivial. 

In the fifth column we indicate whether the different subcycles corresponds to zeros, poles or 
regular points. Many entries required explicit calculations. 

As a special case we can state the following 
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Theorem 3 The sets of zeros and poles ofR[z) are dense on the unit circle. 

To show that the zeros (poles) are dense it suffice to consider points z = exp(2-7rip / q) where q 
is a power of two (three), q — 2 m (q — 3 m ). The function S(z) can thus not be continued outside 
the unit circle. 

We will only speculate what happens for irrational x, that is points that are not periodic 
points of the binary shift map. We know that this map is ergodic and the invariant density is 
uniform. Therefore, the fact the the following integral is zero 

/ log(l - e 2 " x )dx = , (17) 
Jo 

leads us to believe that irrational x corresponds to regular points of S(z). 



4 Concluding remarks 

As we remarked before, our findings concerning the topological zeta functions can be carried over 
directly to general zeta function 

i/CG9,z) = n( 1 -pgr=i9) > ( 18 ) 

only for piecewise linear maps when < /3 < 1 . However for expanding maps, the topological 
zeta function will, after proper rescaling of z, a la eq. (Q), approximate the exact expansion of the 
zeta function and correspond to the fundamental part in a cycle expansion |l9| . It is not clear in 
what sense the curvature corrections of jl9| are small in this case but it is natural to expect the 
natural boundary to prevail. However, the methods used in this paper do not immediately apply 
to the non-linear case so we leave that as a speculation. 

The reasoning above suggests that one could expect a boundary along 

\z\ = {\ff- ) (19) 

when there is no finite Markov partition. (.) is some suitable defined average, presumably geo- 
metric. 

In ref. |2l]| it is shown that for piecewise monotone maps l/((z) is holomorphic inside 

n 

z< lim S np(n\f(xi)\)^ (20) 

n — >oo 

where Xi + \ = f(xi) and the supremum is taken over the initial point x\. The only requirement is 
that the weight, in this case 1/|/'| 1- ' 3 is of finite variation. (For more recent development on this 
subject see ||^, |2^, |24}|.) If there exists a finite Markov partition we know that this holomorphic 
disk may be limited by a pole of 1/C,(z). What we suggest is that in the generic situation with no 
finite Markov partition this guaranteed holomorphic disk has encountered a natural boundary. 

The convergence of cycle expansion depends on the distance between the zero to be determined 
and the nearest singularity of 1/Q(z), zeros close to a natural boundary will be difficult to locate 
accurately. When looking for zeros of the truncated power series expansion of l/((z) one will 
typically find that zeros accumulate along the boundary as the truncation length is increased 



25 . The situation would hardly be improved by applying some continued fraction expansion 



of the zeta function. Such an approach could be a good idea if the (inverse) zeta function is 
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meromorphic, we have encountered such a situation for the preperiodic case. But the presence of 
a natural boundary would presumably lead to accumulation of zeros and poles along the boundary 
of the successive rational approximations of the zeta function. 

I would like to thank Vivianc Baladi and Predrag Cvitanovic for discussions and remarks. 
This work was supported by the Swedish Natural Science Research Council (NFR) under contract 
no. F-FU 06420-303. 



Appendix A 

In this appendix we will prove Theorem 1, that is, we will show that the expansion coefficients 
{bi} in 

n-l 

vttop = - z p) = a - z ) ■ e b ^ ( 2i ) 

p i=0 

are given by 

60 = 1 . (22) 

bi = 6,-ifl, = YY J= i a i 1 < i < n — 1 , 

where the coefficients {ai} are obtained from the finite kneading sequence I x = PC, of length 
n, where P = S1S2 ■ ■ ■ s n _i, according to the rule that ai = 1 if Si = L and a, = —1 if Sj = R. 
The result follows directly from the results of Milnor and Thurston ]20|j . We first recall some 
expressions from [pp| . Let J be an subinterval of the unit interval, and /" the n'th iterate of the 
map. Then one defines 

f 1 f m {x) € J and > 

9(x,JJ m (x))= I -l fm {x) e j and dQsl < o . (23) 

I f m (x)?J 

We use the convention that f°(x) = x. Let the unit interval be divided into I\ and I2 by the 
critical point c (corresponding to the symbol C). For a map with one turning point, the kneading 
matrix is a 1 x 2 matrix (Ki, K2) where 

00 

K i =J2(6(c-,Ii,n-0(c + ,IiJ m ))z m ■ (24) 

m=0 



Perturbation of the critical point as initial value leads to the sequence C ± PC + P of P contains 

■ an odd. We also note that if 1 < m < n — 1 

6m-i (25) 



an even number of R's and C ± PC P if it contains an odd. We also note that if 1 < m < n — 1 

df m (c 



dx 
df m (c-) 
dx 



+b m -i ■ (26) 
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A short calculation now yields 

ifr + if 2 = m - X = 2z ^ m -° m . (27) 

1 — Z n 1 — 2" 

The matrix elements ifi and K 2 are related by 

(l-z)K 1 + (l + z)K 2 = . (28) 
We can now express the kneading determinant 

A = ( M ^ = ESW! (29) 

V 7 1 + 2 2z(l + z) l-z n y ' 

This determinant is related to the inverse topological zeta function 1 / £ top according to |po|| 

n-l 

1/Ctop = D(Z){1 - Z n )(l - Z) = (1 - 2) J2 b mZ m (30) 

m=0 

as claimed. 

If the kneading sequence is infinite (i.e. preperiodic or aperiodic) then the factor (1 — z n ) 
disappears from eq. (|9|) (obvious) and from eq. j30| ) (according to ref. (2Cj) and the only 
alteration of the final result is that n is replaced by oo. 



Appendix B 

Consider the function 

oo 

Z(z = r e 2mx °) = Y[(l - e 2 ™") . (31) 
n=0 

We are interested in the behaviour of 5 when ro — > 1~ . xq is kept fixed and choosen so that the 
sequence {x n } is periodic: x n + m = x n - We rewrite the previous equation as 

oo 

Z(z = r e 2 ™°) = Hf(rf k ) , (32) 

fe=0 

where 

m— 1 

f(r)= l[(l-r 2l e 2 ™") . (33) 
z=o 

First we consider the case |/(1)| > 1, and we will show that |3| — > oo as ro — > 1~. Obviously 
f(r) is an entire analytic function and in particular a continous function so that there exist some 
f such that \f(r)\ > C > 1 for all 1 > r > f. We split up the product 

k— 1 oo k — 1 

S(r e 2 ™«) = J] f(rf k ) • II /(O = ^ ' II /(O ■ (34) 

fc=0 fc= fe fe=0 

We choose r\ = r so that -D(r) is a constant 

oo oo 

w=n/(0=n/(0 > (35) 
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depending only on our choice of f. It is finite because r < 1 and S(z) converges inside the unit 
circle. We now have 

m>D(f)-C k , (36) 

which can be made arbitrary large by choosing k sufficiently large and thus r sufficiently close to 
unity. 

The case when |/(1)| < 1 is worked out in complete analogy it is is shown that |E| — > as 
r ->l- 

The case f(l) — 1 requires just a little more work. First we rewrite eq. ( |34| ) 

k — l k 

Z(roe 2 ™° ) - D(r) ■ JT /(f 2 "*^) = D(f) ■ f[ /(f 2 ^") . (37) 

k=0 k=l 

We want to show that the product converges when k — > oo, and thus ro — » 1 _ . We appeal again 
to the analyticity of f(r). Since /(l) = 1 there exist a constant _E such that 

\f(r) - 1| < £7(1 - r) f < r < 1 . (38) 

The product ( |37| ) converges if 

oo 

1|< oo • (39) 

fc=0 

Using ( pl| ) we get the condition 

f^lKf^^-lK^f^a-f 2 ^) , (40) 

fc=0 fe=0 

which converges since 

(1 -r 2 ~ m ") = 0(-logf 2-" lfc ) fc->oo . (41) 
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